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§ 1. Introduction 

The study of light-matter interactions in dielectric structures with periodic modu- 
lation of the refractive index has developed into a vast research area. At the heart 
of this area is the interplay between Bragg reflections, which block the propagation 
of light in spectral bands known as photonic band gaps (PBGs), and the dynam- 
ical modifications of these reflections by nonlinear light-matter interactions (see 
bibliography compiled by Dowling and Evcritt [2000]). Three- or two-dimensional 
(3D or 2D) PBGs are needed in order to extinguish spontaneous emission in all 
possible directions of propagation, which requires the nontrivial fabrication of 3D- 
or 2D-periodic photonic crystals (Yablonovitch [1987], Yablonovitch [1993]). For 
controlling strictly unidirectional propagation, it is sufficient to resort to PBGs in 
one-dimensional (ID) periodic structures (Bragg refiectors or dielectric multi-layer 
mirrors). Illumination of the periodic dielectric structure at a PBG frequency in 
the limit of vanishing nonlinearity leads to exponential decay of the incident field 
amplitude with penetration depth, at the expense of exponential growth of the 
back-scattered (Bragg-reflected) amplitude. However, this reflection may weaken 
or cease altogether, rendering the structure transparent, when the illumination in- 
tensity and the resulting nonlinearity modify the refractive index so as to shift (or 
even close down) the PBG. The pulsed mode of propagation in nonlinear periodic 
structures exhibits a variety of fundamentally unique and technologically interest- 
ing regimes: nonlinear filtering, switching, and distributed-feedback amplification 
(Scalora, Dowling, Bowden and Bloemer [1994a], Scalora, Dowling, Bowden and 
Bloemer [1994b]). Among these regimes, we have chosen here to concentrate on 
the intriguing solitary waves existing in PBGs, known as gap solitons (GS), and 
solitons propagating near PBGs. 

A GS is usually understood as a self-localized moving or standing (quiescent) 
bright region, where light is confined by Bragg reflections against a dark back- 
ground. The soliton spectrum is tuned away from the Bragg resonance by the 
nonlinearity at sufficiently high field intensities. There is also considerable phys- 
ical interest in finding in the vicinity of a PBG a dark soliton (DS), i.e., a 'hole' 
of a fixed shape in a continuous-wave (cw) background field of constant intensity 
(Kivshar and Luther-Davies [1998]). 

The first type of GS was predicted to exist in a Bragg grating filled with a Kerr 
medium, whose nonlinearity is cubic (Christodoulides and Joseph [1989], Aceves 
and Wabnitz [1989], Feng and Kneubuhl [1993]). Detailed theoretical studies of 
these Bragg-grating (Kerr-nonlinear) solitons (see de Stcrkc and Sipc [1994] for a re- 
view) were followed by their experimental observation (Eggleton, Slusher, de Sterke, 
Krug and Sipe [1996]) in a short (< 10 cm) piece of an optical fiber with a resonant 
Bragg grating written on it. In theoretical considerations of solitons in a Bragg 
grating combined with Kerr nonlinearity, a formidable problem is their stability. 
Current experiments are conducted in short fiber pieces, which do not allow us 
to test the solitons' stability. Approximate and rigorous treatments of the stabil- 
ity problem (Malomed and Tasgal [1994], Barashenkov, Pelinovsky and Zemlyanaya 
[1998], de Rossi, Conti and Trillo [1998]) have demonstrated that the Kerr-nonlinear 
Bragg-grating solitons have instability regions in their frequency-amplitude para- 
metric plane, far from the PBG-center, small-amplitude limit. The generation of a 
slow gap soliton in a Kerr - nonlinear periodic grating via Raman transfer of en- 
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ergy from the pump pulse to the Bragg - resonant Stokes components was recently 
proposed by Winful and Perlin [2000]. 

A Bragg grating with quadratic or second-harmonic-generating (SHG) nonlinear- 
ity can also give rise to a rich spectrum of solitons (Peschel, Peschel, Lederer and 
Malomed [1997], He and Drummond [1997], Conti, Trillo and Assanto [1997]). This 
model has been shown to possess a number of remarkable features. In particular 
(Peschel, Peschel, Lederer and Malomed [1997]), contrary to Kerr-nonlinear Bragg 
gratings, in the case of SHG nonlinearity the gap in which solitons may exist is 
partly empty, only part of it being filled with actually existing soliton solutions. 

While the above temporal-domain models pertain to light propagation in fibers, 
distributions of stationary fields in a planar (2D) nonlinear optical waveguide are 
governed by spatial-domain equations (Stegeman and Segev [1999]). In the planar 
waveguide, a Bragg grating can be realized as a system of parallel 'scratches'. The 
soliton spectrum of this model contains (Mak, Malomed and Chu [1998b]) not only 
the fundamental single-humped solitons but also their two-humped bound states, 
which, quite unusually, turn out to be dynamically stable. Moreover, it possesses 
(Champneys and Malomed [2000]), besides the conventional GSs, also embedded 
solitons, which are isolated solitary-wave solutions within the continuous spectrum, 
rather than inside the gap. 

A principally different mechanism of GS formation has been discovered in a 
periodic array of thin layers of resonant two-level systems (TLS) separated by half- 
wavelength nonabsorbing dielectric layers, i.e., a resonantly absorbing Bragg reflec- 
tor (RABR) (Kozhekin and Kurizki [1995], Kozhekin, Kurizki and Malomed [1998], 
Opatrny, Malomed and Kurizki [1999]). The RABR has been found to allow, for 
any Bragg reflectivity, a vast family of stable solitons, both standing and moving 
(Kozhekin and Kurizki [1995], Kozhekin, Kurizki and Malomed [1998]). As opposed 
to the 27r-solitons arising in self-induced transparency, i.e., resonant field - TLS in- 
teraction in uniform media (McCall and Hahn [1969], McCall and Hahn [1970]), gap 
solitons in RABR may have an arbitrary pulse area (Kozhekin and Kurizki [1995], 
Kozhekin, Kurizki and Malomed [1998]). As shown below, GS solutions have been 
consistently obtained only in a RABR with thin active TLS layers. By contrast, 
an attempt (Akozbek and John [1998]) to obtain such solutions in a periodic struc- 
ture uniformly filled with active TLS is physically unfounded, and fails for many 
parameter values, so that such generalization remains an open problem (see Sec. 
4.3 below). 

An unexpected property of the RABR is that, alongside the stable bright-soliton 
solutions (Kozhekin and Kurizki [1995], Kozhekin, Kurizki and Malomed [1998]), 
this system gives rise to a family of dark solitons (DSs), a large part of which are 
stable (Opatrny, Malomed and Kurizki [1999]). While the existence of stable bright- 
soliton solutions along with unstable DSs is a known feature of uniform SHG media 
(Kivshar and Luther-Davies [1998]), a RABR with thin active layers provides, to 
the best of our knowledge, the first example of a nonlinear optical medium in which 
stable bright and dark solitons exist for the same values of the model's parameters 
(albeit at different frequencies). It is also the first example of the existence of stable 
bright solitons alongside stable cw (background) solutions. 

Potential applications of GSs are based on the system's ability to filter out (by 
means of Bragg reflections) all pulses except for those satisfying the GS dispersion 
condition, as well as to control the pulses' shape and velocity. It would be clearly 
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desirable to supplement these advantageous properties by immunity to transverse 
diffraction of the pulse, i.e., to achieve simultaneous transverse and longitudinal 
self-localization of light in a RABR. This motivates a quest for multi-dimensional 
solitons that are localized in both space and time. 

The concept of optical multi-dimensional spatio-temporal solitons, alias 'light 
bullets' (LBs), was pioneered by Silberberg (Silberberg [1990]), and has since been 
investigated in various nonlinear optical media, with particular emphasis on the 
stability of LBs. For a SHG medium, the existence of stable two- and three- 
dimensional (2D and 3D) solitons was predicted as early as in 1981 (Kanashov and 
Rubenchik [1981]), followed by detailed studies of their propagation and stability 
against collapse (Hayata and Koshiba [1993], Malomed, Drummond, He, Berntson, 
Anderson and Lisak [1997], Mihalache, Mazilu, Malomed and Torner [1998], He 
and Drummond [1998]), Recently, the first experimental observation of a quasi-2D 
spatiotemporal soliton in a 3D SHG sample was reported (Liu, Qian and Wise 
[1999]). The concept of dark LBs was proposed by Chen and Atai [1995]). Stable 
antidark LBs, i.e., those supported by a finite cw background, were predicted in 
a generalized nonlinear Schrodinger equation which contains third-order temporal 
dispersion (Prantzeskakis, Hizanidis, Malomed and Polymilis [1998]). 

As early as in 1984, simulations indicated that self-focusing of spatiotemporal 
pulses in a SIT medium could result in the formation of a quasi-stable vibrating 
object (Drummond [1984]), which was a hint at the possible existence of LBs. It 
has analytically been demonstrated (Blaauboer, Malomed and Kurizki [2000]) that 
uniform 2D and 3D SIT media can indeed carry stable LBs. This investigation 
has been extended to the case of RABR, wherein stable, transversely localized 
SIT solutions combining LB and GS properties are predicted (Blaauboer, Kurizki 
and Malomed [2000]). RABR with any Bragg reflectivity can support stable LBs, 
which are closely related to those in uniform SIT media (Blaauboer, Malomed and 
Kurizki [2000]). 2D LBs supported by a combination of a Bragg reflector with SHG 
nonlinearity were theoretically investigated still earlier (He and Drummond [1998]). 

Finally, we briefly mention a topic that is outside the scope of this review, namely, 
quantum bright solitons, which have been a subject of extensive studies in recent 
years. It has been established that solitons are superpositions of quantum states 
that correspond to clusters of photons bound together. Most of the initial activity 
in this area was devoted to Kerr-nonlinear fiber solitons (Lai and Haus [1989a], Lai 
and Haus [1989b], Wright [1991], Yudson [1985], Kartner and Haus [1993], Cheng 
[1991]). It was shown that optical fiber solitons have quantum analogs which are 
described by the quantum nonlinear Schrodinger equation (QNLSE). The motiva- 
tion was to gain understanding of squeezing effects in soliton propagation (Carter, 
Drummond, Reid and Shelby [1987], Watanabe, Nakano, Honold and Yamamoto 
[1989], Drummond, Carter and Shelby [1989], Haus and Lai [1990], Rosenbluh 
and Shelby [1991] see a review by Sizmann and Leuchs [1999]), two-photon bind- 
ing effects (Deutsch, Chiao and Garrison [1992], Deutsch, Chiao and J. C. Garrison 
[1993]), and fundamental limits imposed on communication systems employing soli- 
tons (Drummond, Shelby, Friberg and Yamamoto [1993], Chiao. Deutsch, Garrison 
and Wright [1993]). Multidimensional quantum solitons were predicted in Kerr and 
SHG waveguides (Drummond and He [1997], Khenmtsyan and Drummond [1998b], 
Kheruntsyan and Drummond [1998a], Kheruntsyan and Drummond [2000]). 

The quantization of GSs has attracted considerable attention as well. A Bethe- 
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ansatz solution was given (Cheng and Kurizki [1995]) for quantum GSs consisting 
of pairwise interacting massive photons (in the effective-mass regime of PBGs prop- 
agation) in Kerr-nonlinear ID periodic grating. A mechanism has been found for 
the creation of two-photon bound states by photons resonantly interacting with 
identical two-level atoms near PBGs in ID periodic structures (Kurizki, Kofman, 
Kozhekin and Cheng [1996]). A Bethe ansatz solution for photons in a PBG ma- 
terial interacting with a single atom was obtained by Rupasov and Singh [1996a], 
Rupasov and Singh [1996b] . It has later been generalized to an extended many-atom 
periodic system, where quantum GSs involving pairs of photons and propagating 
inside a PBG have been found (John and Rupasov [1999]). This subject is of 
potential interest for quantum communications via entangled two-photon states. 

This review starts with a brief survey of solitons in Bragg gratings with cubic and 
quadratic nonlinearities (Ch. 2) and of self-induced transparency (SIT) in uniform 
media and thin films (Ch. 3). It then continues with the derivation of the model for 
phenomena similar to SIT in resonantly absorbing Bragg reflectors (RABR) (Ch. 
4). Bright and dark solitons in RABR are discussed in Ch. 5 and 6, consecutively. 
Light bullets in periodic resonantly-absorbing media are treated in Ch. 7. Finally, 
the prospects for experimental progress in this area are summarized in Ch. 8. 

§ 2. Solitons in Bragg gratings with cubic and quadratic nonlinearities 
2.1. KERR NONLINEARITY 

Because the Bragg grating gives rise to very strong effective dispersion, its combina- 
tion with various optical nonlinearities can create a rich variety of solitons, which, 
in most cases, are gap solitons (GSs), as their intrinsic frequency must belong to a 
gap in the spectrum of linear waves in the Bragg grating. Initially, a Bragg grating 
filled with a Kerr nonlinear medium, whose nonlinearity is cubic, was considered 
by Christodoulides and Joseph [1989], Aceves and Wabnitz [1989], and Feng and 
Kneubuhl [1993]. The corresponding system of two propagation equations for the 
right (forward)- and left (backward)-travcling field envelopes £f{Cj t') and £b{Cj ''") 
with the self- and cross-phase-modulation cubic terms is a generalization of the 
known Thirring model (in the Thirring model proper, the self-phase modulation 
terms are omitted): 



Here the dimensionlcss time r and length C arc related to the physical time t and 
length z as t — At/ (aiUc), and z — 4cC/(aiWc'T-o)j where no is the mean linear 
index of refraction, Ui is its modulation (see Eq. (4.1) below), and u)c is the cen- 
tral frequency of the band gap. Actual values of the electric fields can be obtained 
from the dimensionless quantities £f,b upon multiplication by Wcx/ aino / {ASnx^^'' ) , 
where x^^^ is the third order susceptibility. In Ch. 4.1 we present a general deriva- 
tion of the Maxwell equations in Bragg gratings. 

The existence of solitons in a given model is usually closely related to the modula- 
tional instability of a continuous- wave (cw) solution in the model (Agrawal [1995]): 




(2.1b) 



(2.1a) 
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Figure 1: Experimental evidence for the existence of solitons in a nonlinear optical 
fiber with Bragg grating: the broad pulse is far from a Bragg resonance and has 
passed the fiber without interacting with the grating; the narrow one is a soliton 
shaped by Bragg-resonant interaction with the grating (see further details in the 
work by Eggleton, Slusher, de Sterke, Krug and Sipe [1996]). 

an array (chain) of bright solitons can be generated from the cw solution by modula- 
tional instability. In an optical fiber combining a Bragg grating and Kerr nonlinear- 
ity, all the cw states are modulationally unstable, as observed in a direct experiment 
(which also involved the polarization of light) by Slusher, Spaltcr, Eggleton, Pereira 
and Sipe [2000]. However, as was shown by Opatrny, Malomed and Kurizki [1999] 
(see Ch. 6 below) a RABR gives rise to stable bright solitons coexisting with stable 
cw states, which is a very unusual property. 

Unlike the Thirring model proper, its optical version based on Eqs. (2.1a) and 
(2.1b) is not integrable. Moreover, the equations lack any invariance with respect to 
the reference frame (i.e., the model is neither Galilean nor Lorentz invariant). Nev- 
ertheless, exact single-soliton solutions, which contain two independent parameters, 
viz., an intrinsic frequency ui and the soliton's velocity v, were found in the above- 
mentioned works by Christodoulides and Joseph [1989] and Acevcs and Wabnitz 
[1989], following the pattern of the Thirring-model solitons. At v = 0, the quies- 
cent (standing) solitons completely fill the gap — 1 < w < +1 in the spectrum of the 
linearized equations (2.1a) and (2.1b), and the velocity takes all the possible values 
— 1 < t; < +1 (recall that, in the present notation, the maximum group velocity of 
light is 1). Solitons with small frequencies have small amplitudes and are close to 
the classical nonlinear Schrodinger (NLS) solitons, while the ones with values of ui"^ 
closer to 1 are strongly different from their NLS counterparts; in particular, they 
are chirped, i.e., they have a nontrivial intrinsic phase structure. 

Many properties of these Bragg-grating solitons, as they are frequently called, 
were reviewed by de Sterke and Sipe [1994] (see also the special issue of Optics 
Express edited by Brown and Eggleton [1998]). Detailed theoretical studies were 
finally followed by their experimental observation by Eggleton, Slusher, de Sterke, 
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Knig and Sipe [1996] in a short (< 10 cm) piece of an optical fiber with a resonant 
Bragg grating written on it. This experiment requires the use of very powerful light 
beams, with an intensity comparable to the fiber's breakdown threshold. Despite 
this difllculty, propagating pulses with soliton-like shapes have been detected in the 
experiment, see Fig. 1. 



It should be stressed that the above-mentioned exact soliton solutions to Eqs. 
(2.1a) and (2.1b) take the simplest form in the case v = 0, corresponding to a pulse 
of standing light. In reality, the soliton observed by Eggleton, Slusher, de Sterke, 
Krug and Sipe [1996] was moving at a considerable velocity. Generation and detec- 
tion of zero-velocity solitons in nonlinear Bragg gratings remains an experimental 
challenge. 



Theoretical considerations of solitons in Bragg gratings with Kerr nonlinearity 
face the tough problem of soliton stability. The first approach to the stability prob- 
lem, developed by Malomed and Tasgal [1994] was based on the variational ap- 
proximation. Although this approach is not rigorous, it clearly demonstrates that, 
sufficiently far from the above-mentioned NLS (small-frequency small-amplitude) 
limit, Bragg-grating solitons can be unstable. Later, a rigorous analysis of the same 
stability problem was developed by Barashenkov, Pelinovsky and Zemlyanaya [1998] 
and by de Rossi, Conti and Trillo [1998]. It has been demonstrated that the Bragg- 
grating solitons indeed have instability regions in their parametric (w, v) plane, 
which are close to those originally predicted by means of the variational approxi- 
mation. 



The derivation of the standard equations (2.1a) and (2.1b) from the imderlying 
Maxwell's equations (see Ch.. 4.1) neglects all the linear terms with second-order 
derivatives (which account for the material dispersion and/or diffraction in the 
medium). As shown by Champncys, Malomed and Friedman [1998], taking the 
second derivatives into account drastically changes the soliton content of the model 
(although the propagation distance necessary to observe the change of the solitons' 
shape may be much larger than that available in current experiments). The rel- 
atively simple soliton solutions generated by Eqs. (2.1a) and (2.1b) immediately 
disappear after the inclusion of the second-derivative terms; instead, three com- 
pletely new branches of soliton solutions emerge, provided that the coefficient in 
front of the second-derivative terms exceeds a certain minimum value. Very re- 
cently, it has been shown by SchoUmann and Mayer [2000] that these new branches 
are, generally, also subject to an instability. 



Finally, we mention a dual-core nonlinear optical fiber, each core containing 
a Bragg grating, which is described by four equations (Mak, Malomed and Chu 
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[1998a]), 



.dSpi .dSpi f^ic i2 \c i2 



+ Sbi + \£f2 - 0, (2.2a) 
+ Efi + \£bi = 0, (2.2b) 

.dEF2 , .QEfi , /l|c. |2 , ic |2 1 c> 



dr dC V2 

+ 5S2 + A£fi = 0, (2.2c) 

.d£B2 .d£B2 , |2 , ic |2 1 c 
+ T: ^B2 + ^F2 Kb2 



ar aC V2 

+ £fi + \£bi = Q, (2.2d) 

where the subscripts 1 and 2 stand for the core number, the dimensionless quantities 
T, C, £b,F;1,2 have the same meaning as in Eq. (2.1), and A is the real hnear coupUng 
between the cores. Due to this Hnear coupling, a soliton in this system necessarily 
has its components in both cores. As shown by Mak, Malomed and Chu [1998a], 
an obviously symmetric soliton with equal components in the two cores, £fi = £f2 
and £bi = £b2, becomes unstable when its energy exceeds a certain threshold. 
The instability gives rise to a pitchfork bifurcation that generates a pair of stable 
asymmetric solitons which are mirror images of each other. 

2.2. QUADRATIC NONLINEARITIES 

A Bragg grating with quadratic (sccond-harmonic-gcncrating, SHG) nonlinearity 
gives rise to a vast variety of solitons. This model was introduced independently and 
simultaneously by Pcschcl, Pcschcl, Ledercr and Malomed [1997], Ho and Drum- 
mond [1997], and Conti, Trillo and Assanto [1997]. The model includes four fields, 
namely, the fundamental- and second-harmonic components of the forward- (F-) 
and backward- (B-) traveling waves. In a normalized notation, the corresponding 
system takes the form 

i^+i^+x£F + £*FSF+£B = 0, (2.3a) 



dr dC 
.IBQf .OGf 
V OT oQ 
. d£B _ . d£B 

' dr ' dc 

.IdGB .OGb 



\-X~^Qf+£% + xGb = 0, (2.3b) 
+ x£b+£*bOb+£f = 0, (2.3c) 

. ^ ^X~^Gb+£1 + ^Gf = 0, (2.3d) 

where £ and G are the fundamental- and second-harmonic fields, t; > is the 
group velocity of the second harmonic relative to the fundamental harmonic, x is 
a (generally complex) coupling F-B coefficient in the second harmonic, given that 
the coupling constant at the fundamental harmonic is normalized to 1, and the real 
parameter x determines the phase mismatch of the two harmonics. 
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Figure 2; An example of a dynamically stable double-humped soliton in the model 
based on Eqs. (2.3a) - (2.3d), combining a Bragg grating and second-harmonic- 
generation nonlinearity. Shown in this figure is the intensity of the forward 
fundamental-harmonic component I^fI^ (see Peschel, Peschel, Lederer and Mal- 
omed [1997]). 

Analytical and numerical consideration of solitons in this model has revealed a 
number of nontrivial features. In particular, Peschel, Peschel, Lederer and Malomed 
[1997] have found that, contrary to the model (2.1) with Kerr nonlinearity, in the 
system (2.3) the gap in which solitons may exist has voids, so that only a part of 
it is filled with actually existing soliton solutions. Another noteworthy finding of 
the same work is that not only single-humped fundamental solitons, but also their 
two-humped bound states, which are normally unstable in other versions of the 
SHG models, arc dynamically stable in the four- wave model (2.3). An example of a 
two-humped soliton which turns out to be fairly stable in simulations of the system 
(2.3) is shown in Fig. 2. A detailed description of gap solitons in the four- wave 
model combining a Bragg grating and SHG nonlinearity can be found in the review 
by Etrich, Lederer, Malomed, Peschel and Peschel [2000] on solitons in various SHG 
media. 

A different model based on resonant Bragg reflection and SHG nonlinearity can 
be formulated in the spatial domain [unlike the time-domain models (2.1) and 
(2.3)], for stationary fields in a planar (2D) nonlinear optical waveguide. In such a 
waveguide, a ID Bragg grating is realized as a system of parallel scores ('scratches'). 
The simplest model of this type involves three waves (Mak, Malomed and Chu 
[1998b]): 

1^+1^+82 + 838^ = 0, (2.4a) 



i^-i^+S,+S-iSl = 0, (2.4b) 
2i^-q£s + D^+£iS2 = 0. (2.4c) 



Here C, and x are the propagation and transverse coordinates, respectively; the fields 
£1 and £2 are two components of the fundamental harmonic that are transformed 
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into each other by resonant reflections on the ID Bragg grating, £3 is the second- 
harmonic component, and D is an effective diffraction coefficient for the second 
harmonic. The wave vectors ki 2,3 of the three waves are related by the resonance 
condition, ki+k2= ks, the real parameter q accounts for a residual phase- mismatch. 
The configuration corresponding to this model assumes that the second harmonic 
propagates parallel to the Bragg grating (which has the form of the above-mentioned 
scores). It is therefore necessary to take into account the diffraction of this com- 
ponent, while for the two fundamental harmonics the effective diffraction induced 
by resonant Bragg scattering is much stronger than normal diffraction, which is 
neglected here. 

The soliton spectrum of this model is fairly rich. It contains (Mak, Malomed 
and Chu [1998b]) not only fundamental single- humped solitons but also their two- 
humped boimd states, some of which, as in the case of the four- wave model (2.3), 
may be dynamically stable. A rigorous stability analysis for various solitons in the 
model (2.4), based on computation of eigenvalues of the corresponding linearized 
equations, was performed by Schollmann and Mayer [2000] . This analysis has shown 
that some of these solitons, although quite stable in direct dynamical simulations, 
are subject to a very weak oscillatory instability, whereas other solitons in this 
model are stable in the rigorous sense. 

The three- wave model (2.4) possesses (Champneys and Malomed [2000]), besides 
the traditional GSs, numerous branches of embedded solitons: isolated solitary- 
wave sohitions existing within the continuous spectrum, rather than inside the gap. 
Solutions of this kind appear also when the second-derivative terms are added to 
the generalized Thirring model (2.1). 

Finally, the four-wave model (2.3) with quadratic nonlinearity can be extended to 
the two- and three-dimensional cases, by adding transverse diffraction terms to each 
equation of the system. Physically, this generalization corresponds to spatiotempo- 
ral evolution of the fields in a two- or three-dimensional layered medium. Because, 
as is well known, quadratic nonlinearity does not give rise to wave collapse in any 
number of physical dimensions, the latter model can support stable spatiotempo- 
ral solitons, frequently called light bullets. Direct numerical simulations reported 
by He and Drummond [1998] have confirmed the existence of stable 'bullets' in a 
multidimensional SHG medium embedded in a Bragg grating. 

§ 3. Self-Induced Transparency (SIT) in Uniform Media and Thin 
Films 

3.1. SIT IN UNIFORM MEDIA 

Self-induced transparency (SIT) is the solitary propagation of electromagnetic (EM) 
pulses in near-resonant atomic media, irrespective of the carrier-frequency detuning 
from resonance. This striking effect, which is of paramount importance in nonlinear 
optics, was discovered by McCall and Hahn [1969], McCall and Hahn [1970]. If the 
pulse duration is much shorter than the transition (spontaneous-decay) lifetime 
(Ti) and dephasing time (T2), then the leading edge of the pulse is absorbed, 
inverting the atomic population, while the remainder of the pulse causes atoms to 
emit stimulated light and thus return the energy to the field. When conditions for 
the process are met, it is found that a steady-state pulse envelope is established 
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and then propagates without attenuation at a velocity that may be considerably 
less than the phase velocity of light in the medium. 

We start with the Hamiltonian for a single atom in the field, 

H=^u^.E-d, (3.1) 

where 

w=\e){e\-\g){g\, (3.2) 

is the atomic inversion operator, wq is the atomic transition frequency, \g) and |e) 
denote the atomic ground and excited states, respectively, E is the electric field 
vector and d is the atomic dipole-moment operator. We take the projection on the 
field direction, so that E • d = Ed, where 

d^|(p + pt), (3.3) 

fj, being the dipole moment matrix element (chosen real) and 

P^2|5)(e| (3.4) 

atomic polarization operator. 
We express the electric field at a given point by means of the Rabi frequency Q 

as 

E = —{ne-^"^ + fi*e^'^=*). (3.5) 



for the atomic po- 



The Heisenberg equations of motion dA/dt = l/{ih) 
larization and inversion operators (3. 4), (3. 2), yield the Bloch equations for their 
expectation values (c-numbers) P and w, respectively 

dtP{z,t) = w{z,t)Sl — i{u)o — ijOi2)P, (3.6a) 
dtw{z,t) = ~[P*{z,t)Q.+ c.c.]. (3.6b) 

The Maxwell equations (Newell and Moloney [1992]) reduce in the rotating- wave 
and slow- varying approximations to 



no dz dt 



where 



ro = ^W^, (3.8) 

is the cooperative resonant absorption time, qq being the TLS density (averaged 
over z), and no is the refraction index of the host media. 
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In the simplest case, when the driving field is in resonance with the atomic 
transition, ujq — lo^, the Bloch equations (3.6) can be easily integrated and the 
Maxwell equation (3.7) then reduces to the sine-Gordon equation 

- = -sm.6 (3.9) 



dCdf 

for the 'rotation angle', 



t 

ndt', (3.10) 

oo 



in terms of the dimensionless variables f = {t — noz/c)/To and ( = uoz/cto. 

This sine-Gordon equation is known to have solitary-wave solutions, for which 
the total area under the pulse is conserved and equal to 27r - the so called pulse-area 
theorem by McGall and Hahn [1969], McGall and Hahn [1970]: 

0(C, f ) = (ro)-^Aosech [/3(C - vf)] , (3.11) 

where pulse width f3 is an arbitrary real parameter uniquely defining amplitude 
Aq = 2/(3 and group velocity v = 1/ (fi of the soliton. 

Since its inception, SIT has become an active research area with many practical 
applications, for which we refer readers to excellent reviews by Lamb- Jr. [1971], 
Poluektov, Popov and Roitberg [1975], Maimistov, Basharov and Elyutin [1990] 
and references therein. In this section we will only briefly discuss results which 
are pertinent to the present review, such as SIT in thin films and collisions of 
counter-propagating SIT solitons. 



3.2. SIT IN THIN FILMS 

The interaction of light with a thin film of a nonlinear resonant medium located at 
the interface between two linear media has been described by Rupasov and Yudson 
[1982] and Rupasov and Yudson [1987], who have shown that a nonlinear thin film 
of TLS can be a nearly ideal mirror for weak pulses, but transparent for pulses of 
sufficient intensity. The problem of light pulse transmission through the nonlinear 
medium boundary has been studied under conditions of coherent interaction with 
the matter. The system can be described by a set of nonlinear Maxwell-Bloch-like 
equations which effectively take the presence of the refiected wave into account by 
imposing boundary conditions on the electromagnetic fields at the interface. It has 
been shown (Rupasov and Yudson [1987]) that these equations are exactly inte- 
grable by the inverse scattering method, and 27r-soliton-pulse transmission through 
the film has been studied. 

If the atomic density is such that on average there are more than one atom per 
cubic resonant wavelength, then near-dipole-dipole (NDD) interactions, or local- 
field effects, can no longer be ignored, contrary to the case of more dilute media. 
NDD effects necessitate a correction to the field that couples to an atom in terms 
of the incident field and volume polarization (Bowden, Postan and Inguva [1991], 
Scalora and Bowden [1995]). This effect can give rise to bistable optical transmis- 
sion of ultrashort light pulses through a thin layer consisting of two-level atoms 
(Basharov [1988], Benedict, Malyshev, Trifonov and Zaitsev [1991]): the local-field 
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correction leads to an inversion-dependent resonance frequency, and generates a 
new mechanism of nonlinear transparency. When the excitation frequency is some- 
what larger than the original resonant frequency, the transmission of the layer 
exhibits a transient bistable behavior on the time scale of superradiancc (Basharov 
[1988], Benedict, Malyshev, Trifonov and Zaitsev [1991]). It was shown that if an 
ultrashort pulse is allowed to interact with a thin film of optically dense two-level 
systems, the medium response is characterized by a rapid switching effect (Cren- 
shaw, Scalora and Bowden [1992], Crenshaw and Bowden [1992]). This behavior 
is more remarkable than the response of conventional two-level systems, because 
the medium can only be found in one of two states: cither fully inverted, or in the 
ground state, depending (quasi-periodically) on the ratio between the peak field- 
strength and the NDD coupling strength. This feature was found to be impervious 
to changes in pulse shape, and independent of the pulse area (Crenshaw, Scalora 
and Bowden [1992]). 

Passage of light through a system of two thin TLS films of two-level atoms has 
been considered by Logvin and Samson [1992], Logvin and Loiko [2000] who have 
shown that if the distance between the films is an integer multiple of the wavelength, 
then the system is bistable. Self-pulsations, i.e., periodically generated output, arise 
if an odd number of half- wavelengths can be fitted between the films and absorption 
in the medium is insignificant. In general, the dynamics admits both regular and 
chaotic regimes. 

3.3. COLLISIONS OF COUNTERPROPAGATING SIT SOLITONS 

Situations in which it is necessary to consider the interaction of incident (forward) 
and reflected (backward) light waves include: intrinsic optical bistability (Inguva 
and Bowden [1990]), dynamics of excitations in a cavity (Shaw and Shore [1990]) 
and collisions of counterpropagating SIT solitons (Afanas'ev, Volkov, Dritz and 
Samson [1990], Shaw and Shore [1991]). The field in such problems is represented 
as a superposition of forward- and backward-traveling waves. The; atomic response 
to this field is determined by solving the Bloch equations (3.6) in the rotating- 
wave approximation. The population inversion w{z,t) and polarization P may 
be represented by a quasi-Fourier expansion over a succession of spatial harmonic 
carriers and slow varying envelopes, entangled in a fashion which leads to an infinite 
hierarchy of equations. The truncation of this hierarchy can only be justified by 
phenomenological arguments, such as atom movement in an active atomic gas. 

When the forward (F-) and backward (B-) wave pulses overlap in space and time, 
the resulting interference pattern of nodes modifies the atomic excitation pattern. 
The spatial quasi-Fourier expansion provides an efficient way of treating the spatial 
inhomogeneities of the response in those regions where the F- and B-pulses overlap, 
each successive Rabi cycle increasing the number of terms that contribute to the 
expansion (Shaw and Shore [1990]). 

Collisions of optical solitons produce observable effects on both the atoms and 
the pulses. The overlap of two counterpropagating pulses can produce an appre- 
ciable spatially-localized inversion of the atomic population, thus causing optical 
solitons to lose energy. It was found that, whereas large-energy solitons passed 
freely through each other, solitons whose initial energy fell below a critical value 
were destroyed by collisions. In addition, the residual atomic dipole, created by the 
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excitation, acts as a further source of radiation. This radiation appears as an oscil- 
lating tail on the postcoUisional pulses and, over longer time scales, as fluorescence 
(Afanas'ev, Volkov, Dritz and Samson [1990], Shaw and Shore [1991]). 



§4. SIT in Resonantly Absorbing Bragg Reflectors (RABR): The 
Model 



4.1. MAXWELL EQUATIONS 

Let us assume (Kozhekin and Kurizki [1995], Kozhekin, Kurizki and Malomed 
[1998], Opatrny, Malomed and Kurizki [1999]) a one-dimensional (ID) periodic 

modulation of the linear refractive index n{z) along the z direction of the electro- 
magnetic wave propagation (see Fig. 3). The modulation can be written as the 
Fourier series 



n^iz) = ng[l + oi cos{2kcz) + 02 cos(4fccz) -I- . . • ], 



(4.1) 



where ng, aj and kc arc constants, and the medium is assumed to be infinite and 
homogeneous in the x and y directions. 

The periodic grating gives rise to photonic band gaps (PBGs) in the system's 
linear spectrum, i.e., the medium is totally reflective for waves whose frequency is 
inside the gaps. The central frequency of the fundamental gap is lUc ~ kcc/no, c 
being the vacuum speed of light, and the gap edges are located at the frequencies 



(^1,2 =a;c(l±ai/4), 



(4.2) 



where ai is the modulation depth from Eq. (4.1). 

We further assume that very thin TLS layers (much thinner than l/kc), whose 
resonance frequency ujq is close to the gap center Uc, are placed at the maxima of 
the modulated refraction index. In other words, the thin active layers are placed 
at the points z\a,yei such that cos(A;c^iayer) = ±1. 

We shall study the propagation of electromagnetic waves with frequencies close 
to oJc through the described medium. Let us write the Max;well equation for one 
component of the field vector propagating in the z direction as 



dz^ 



n\z) 



dt^ 



dt^ 



(4.3) 



with the refraction index n modulated as in Eq. (4.1), E being the electric field 
component and P„; the non-linear polarization. We use the substitution 



E = [Spiz, t)e'''-^ + Ssiz, i)e-*'==^] e'^-* + c.c. 



(4.4) 



with ujc satisfying the dispersion relation uquJc = kcC and Sp and £3 denoting the 
forward and backward propagating field components. We work in the slowly- varying 
envelope approximation 



(4.5a) 
(4.5b) 



d^£B,F 

dz"^ 


< 


OZ 


d^£B,F 

dt^ 




d£B,F 

^'^ dt 
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Figure 3: Schematic description of the periodic RABR and of the decomposition 

of the electric field into modes E+ and E_. The shading represents regions with 
different index of refraction; the darker the shading the larger n is. The black 
regions correspond to the TLS layers. The upper solid curve represents the elec- 
tric field, the lower solid curves correspond to the components Re(S+)cosfccZ and 
— Im(S_)sinfcc-2; the dashed curves are the envelopes Re(I]+) and — Im(I]_). The 
vertical dotted lines denote the positions of the TLS. 



Substituting (4.4) into (4.3), using (4.5a) and (4.5b), multiplying by e*(T'=^^+'^<'') 
and averaging over the wavelength A = 27r/fcc and the period T = 2it /uc we get 

+ — = —^^B + T^P-, (4.6a) 



no dz dt 4 
no oz at 4 2/xTo 



where we define 



P± = -^/^e^(^'=-+-*)\ , (4.7) 

with jjL being the dipole moment (Ch.3, and tq a constant chosen here to be the 
medium absorption time [see Eq. (3.8)]. The averaging is defined as 

{...),^,^^ll...dtdz. (4.8) 

We express the field components £f,b by means of the dimensionless quantities S± , 

£:.,. = ^(E,±E_), (4.9) 

so that the electric field E = E{z, t) is 

E{z,t) = hiiiToy^ (Re [j:+{z,t)e-"^'*] coskcZ 

-Im [Y,_{z,t)e-^''*]smkcz) . (4.10) 

To obtain the equations of motion in the most compact form, it is convenient to 
introduce the dimensionless time r, coordinate (, and detuning S as 

T = t/ro, C = (no/cTo) z, S = {u>o - u>c)to, (4-11) 
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where Tq is the characteristic absorption time of the field by the TLS medium as 
defined by Eq. (3.8). Substituting Eq. (4.9) into (4.6a) and (4.6b), using (4.11) and 
(3.8), then difi'erentiating the equations again with respect to ( and r, we arrive, 
after algebra, at the form 



and 



+ ^ (P+ +P-) + §^ {P+ - P-) , (4.12a) 



-^(P+-P_)-^(P++P_). (4.12b) 

Here the dimensionless modulation strength r] is the ratio of the TLS absorption 
distance Zabs = toc/uq to the Bragg reflection distance Zrefl = 4:c/{aiu)cno), which 
can be expressed as 

r] = li,hs/lreR = aiUcTo/4:. (4.13) 

The equations for the electric field components T,± can be solved once we know 
the averaged polarization P±. To find P±, we express the polarization P„j at a 
given point as the dipole moment density 

Pnl = 47Tg{d) ^27rgti(^{P) + {P^)'^ 

= -2TTigfi (Pe-^'^<=* - p*g«-^ct) ^ (4^4) 

where g is the number of the two-level atoms in a unit volume. Neglecting the time 
derivatives of P with respect to those of e~"^<=*, we can write for the polarization 
derivative 

^ = 27riojlgfi (Pe--=* - P*e-=*) , (4.15) 
so that P± on the right-hand side of Eqs. (4.6a), (4.6b) read as 



{gPe^'^"'')^. (4.16) 



To determine the evolution of P± , we have to make certain approximations. Let 
us first consider the situation when the TLS's are confined to infinitely thin layers. 
We will subsequently study the influence of finite width of the layers. 

4.2. TWO-LEVEL SYSTEMS (TLS) IN INFINITELY THIN LAYERS 



Let us now assume that the atomic density g is concentrated in zero-width layers 



located at zj such that 
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i.e., it is described by 



Q=^-^Y.^{z-Zj), (4.18) 

3 

where Qo is the bulk density averaged over the whole wavelength. We also assume 
that the spatial dependence of P is to a good approximation anti-periodic with 
respect to A/2, i.e. P{z+X/2) « — P{x). This is in agreement with the approximate 
anti-periodicity of the electric field with A/2. Denoting by Po the value of P in 
even- numbered layers (the value in the odd- numbered layer being — -Pq)) we get the 
spatial average in (4.16) as 

,2 „ Jl 



Note that P+ = P_ is the consequence of the zero width of the layers. Due to the 
choice of tq as in Eq. (3.8), we obtain the simple relation 

P± = Po. (4.20) 

Using this expression in (4.12a) and (4.12b), we obtain the following evolution 
equations for S-t 

dP 

-- -7]'^E++2ir]P + 2 — , (4.21a) 

OT 















5r2 





Qp 

-- -ri'^_-2—, (4.21b) 

where we have omitted the index of Pq, for simplicity. 

The equations for the polarization P and inversion w in the even- numbered layers 
can be obtained from Eqs. (3.6) by substituting for f2 (combining Eqs. (4.10) and 
(3.5)) 

VI = To"^ (i;+ cos kcz + sin kcz) (4.22) 

and applying Eq. (4.17) at the positions of these layers. Expressing the detuning 
as in Eq. (4.11), we obtain the equations 

dP 

— = -i5P + J:+w, (4.23a) 

OT 

^ = -Re (S+P*). (4.23b) 

The set of equations (4.21) and (4.23) was first derived by Mantsyzov and 
Kuz'min [1984] and Mantsyzov and Kuz'min [1986] for the particular case of a 
periodic array of thin TLS layers without modulation of the linear index of refrac- 
tion, i.e., r] = 0. In this case, these equations can be reduced to the sine-Gordon 
equation (3.9) for the area of the 'forward' wave 

^{z,t) = I Y,+ {z,t')dt' . (4.24) 
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It follows from (4.24) that the 27r solitary wave (SIT), associated with the sum of 
the forward and backward propagating waves, is an exact solution to the coupled- 
mode Maxwell-Bloch equations for a periodic stack of resonant thin films satisfying 
the Bragg condition. These solitary wave solutions are referred to as "two-wave 
solitons" (Mantsyzov and Kuz'min [1986]). 

If the Bragg condition is not exactly satisfied, then the system can exhibit a 
rich, multi-stable behavior, which has been studied numerically (Mantsyzov [1992], 
Lakoba and Mantsyzov [1992], Lakoba [1994]). In the case of a slight violation of the 
Bragg condition at the exact two-level resonance, an analytical solution for the 27r 
gap soliton with a phase modulation had been obtained (Mantsyzov [1995]). These 
studies have also revealed the existence of an oscillating pulse, whose amplitude 
and velocity sign change periodically. 

Combining Eqs. (4.23a) and (4.23b), one can eliminate the TLS population 
inversion: 

w = ±Vl - (4.25) 

Without the field-induced polarization, the TLS population is not inverted (w = 
— 1), hence the lower sign must be chosen in Eq. (4.25). Thus, the remaining 
equations for E+ and P form a closed system, 

rfi:+ + 2i{r] - 6)P 

2^/1 -|P|2 S+, (4.26a) 
iSP - sj\ - |P|2 S+, (4.26b) 

and S_, the field component driven by dP/dC,, can then be found from Eq. (4.21b). 



dP 

a7 



4.3. FINITE WIDTH OF TLS LAYERS 

So far, we have assumed that the TLS layers are infinitesimally thin. We now 
proceed to estimate effects of non-zero width of the active layers, which represent 
more realistic physical situations. We still assume the width of the active layers to 
be small in comparison to the wavelength. This allows us to expand the polarization 
as a Taylor series in the position within the layer, and consider only terms up to 
the second order. Averaging the polarization over the entire wavelength yields the 
source term of the Maxwell equations (4.6). 

Let us assume that the TLS density is given, instead of Eq. (4.18), by 




(4.27) 



where ^iq is the density averaged over the whole wavelength and the function f{z) 
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/(^)>0, f{z + \/2) 

f{z) dz 



L 



Optical Solitons in Periodic Media 
1, 



3 4 



2=3 +T 



■3 4 



z f{z) dz 
z^ f{z) dz 



= 0, 



r 



(4.28) 



where 7 <C 1 is a dimensionless parameter describing the thickness of the layers. 
For a rectangular function / of width D this parameter is 



7 = 



(4.29) 



To calculate (gPe^^'^"^) ^ of (4.16), we expand the spatial dependence of P and 
exp{±ikcz) in Taylor series, 



^2i + i 



A 

22j - 4 



P0 + P^{z-Z2j) + -P^'{z-Z2jf- 
1 ± ikc{z-Z2j) - -^{Z-Z2jf + ■ ■ 



dz, 



(4.30) 



where Po, ^o, and Pq' refer to the values of P and its spatial first and second 
derivatives at the positions of the even- numbered layers. Neglecting higher than 
quadratic terms and using Eqs. (4.28), we arrive at 



° ^ ^ > 2fc2 ^ fc 2 



so that 



P±^Po + Y 



2fc2 k 2 



(4.31) 



(4.32) 



Using these values in (4.16), substituting them into (4.12a) and (4.12b), and defining 

Po = P{Z2,), (4.33) 
. , ,434, 



Pc 



1^-llp 

fc2 dz'^ 



(4.35) 



-Z2j 
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we obtain the field equations of motion in the form 



OP 



dPc BPb 



irjPc + 



dr 



,dPo 



+7^ 



-irjPB - 



dPc dPs 



dr 
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(4.36a) 



(4.36b) 



To obtain the equations of motion for the atomic parameters, we use Eq. (4.22) 
for Q in (3.6) and rewrite the Bloch equations at a point z as 

dP 

— = -i5P+[cos{kcz)'E+ + ism{kcz)'E_]w, (4.37a) 

OT 

d 1 

= -- [cos(A:c^)S+ +isin(A:c^)5]_]P* + C.C.. (4.37b) 

OT 2 

Considering that the spatial derivatives of S± are much smaller than those of 
sin(fcc2;) and cos(fcc-z), and calculating the first and second derivatives at the posi- 
tions of the even-numbered layers, we get 



and 



dPp 
dr 

dn 

dr 
dPl 
dr 



dwp 
dr 

dr 

~d7 



-iSPo + S+wo, 



--S+P*+c.c., 

-^S_Po*-i5]+P-+c.c., 
^E+Po* - ifccS-P^* - ^E+P- + c.c. 



(4.38) 
(4.39) 
(4.40) 

(4.41) 

(4.42) 
(4.43) 



Using the definitions (4.33) - (4.35) and defining 



Wo 
W2 



W{z2j), 

i dw 
k dz 



k^ dz"^ 



(4.44) 
(4.45) 

(4.46) 



Z=Z2j 
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we then arrive at the equations of motion 



-iSPB + 2w;oS- - 2w;iE+, 



-i6Po +wo'S+, 



(4.47) 



(4.48) 



dPc 



-iSPc + t«2S+ + 2wiT, 



(4.49) 



dr 
dwQ 

dr 
dwi 

dW2 

'd7 



--E+Po*+c.c., 
1 1 



-S.P^+c.c. 



c.c. 



(4.50) 



(4.51) 



(4.52) 



Equations (4.36a), (4.36b) and (4.47) - (4.49) with their complex conjugates and 
Eqs. (4.50) - (4.52) form a closed set of 13 equations for the variables S+, S_, Pq, 
Pb, Pc and their complex conjugates and wq, wi and W2, i.e., 13 real variables 
together. The equations are parametrized by 3 real parameters ry, 5 and 7. Note 
that for 7 = Eqs. (4.36a), (4.36b), (4.47) and (4.50) are identical to Eqs. (4.21a), 
(4.21b), (4.23a) and (4.23b), respectively, with Pq standing for P. Even though the 
number of equations and variables has now increased, they are still relatively easy 
to solve numerically, and can realistically express the properties of the system over 
rather long times. 

The case of a uniform active medium embedded in a Bragg grating calls either 
for a solution of the full second- order Maxwell equation, without the spatial slow 
varying approximation, or for a soluton of an infinite set of coupled equations for 
Pi and Wi as in the case discussed in Ch.3.3. This makes the present analysis prin- 
cipally at variance with that of Akozbek and John [1998], where the slow varying 
approximation for atomic inversion and polarization is assiimcd, thus arbitrarily 
truncating the infinite hierarchy of equations to its first two orders. 

4.4. ENERGY DENSITIES 

To reveal the physical meaning of the quantities S± and P, we express the energy 
density of the electromagnetic field as 



that of the TLS excitations (considering the limit of infinitely thin TLS layers) as 



(4.53) 




(4.54) 



and the energy density of the TLS-field interaction as 



Wi = {l/2)hpoT^^ Im (S+P*) . 



(4.55) 



From Eq. (4.53) we conclude that and are proportional to the number 

of photons per TLS (atom), in the standing- wave symmetric and anti-symmetric 
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Figure 4: The RABR dispersion curves (dimensionless frequency x versus dimen- 
sionless wavevector fc) at = 0.5 and S = —0.2. The soHd hues show the dispersion 
branches corresponding to the 'bare' (nonintcracting) grating, while the dashed and 
dash-dotted hues stand for the dispersion branches of the grating 'dressed' by the 
active medium. The frequency bands that support the standing dark and bright 
sohtons arc shaded. The arrow indicates a complete gap, where no field propagation 
of any kind can take place. 



modes, whose anti- nodes and nodes, respectively, coincide with the active layers 
(see Fig. 3). Since the interaction time tq (see Eq. (3.8)) is usually much larger 
than the optical period 2^/(^0, the interaction energy is negligible in comparison 
with the energies of the field and atomic excitations. 

4.5. LINEARIZED SPECTRUM 

To reach general understanding of the dynamics of the model, the first necessary 
step is to consider the spectrum produced by the linearized version of Eqs. (4.21b), 
(4.26a) and (4.26b), which describe weak fields in the limit of infinitely thin TLS 
layers. Setting w = —1, and 

S+ = Ae'^''^-^''\ (4.56a) 
E_ = Be'^'"^-^^\ (4.56b) 
P = Ce'^^^-'^^\ (4.56c) 

we obtain from the linearized equation (4.26b) that C = i{S — x)~^A. Substitut- 
ing this into Eqs. (4.21b) and (4.26a), we arrive at the dispersion relation for the 
wavenumber k and frequency % in the form 

(x' V^)ix -S)x {(x - S) [x^ -K^-{2 + 77^)] +2{fj-d)}= 0. (4.57) 

Different branches of the dispersion relation generated by Eq. (4.57) arc shown 
in Fig. 4. The roots x = + (corresponding to the solid lines in Fig. 4) 

originate from the driven equation (4.21b) and represent the dispersion relation 
of the Bragg reflector with the gap |x| < f] (cf. Eq. (4.2)), that does not feel 
the interaction with the active layers. The degenerate root x = is trivial, as it 
corresponds to the eigenmode (4.56c) with A = B = 0. Important roots are those 



24 



Optical Solitons in Periodic Media 



given by the curly brackets in Eq. (4.57) (shown by the dashed and dash-dotted 
lines in Fig. 3), since they give rise to nontrivial spectral features. They will be 
shown below to correspond to bright or dark solitons in the indicated (shaded) 
bands. 

The frequencies corresponding to fc = are 

Xo = V and xo,± = ~{v-6)± ^27^^^, (4-58) 

while at A:^ oo the asymptotic expressions for different branches of the dispersion 
relation are x = and x = S + 2{ri — S)k~^. Thus, the linearized spectrum 
always splits into f,wo gaps, separated by an allowed band, except for the special 

case, T] = r]o = + \Ji + j^^, when the upper gap closes down. The upper and 
lower band edges are those of the periodic structure, shifted by the induced TLS 
polarization in the limit of a strong reflection. They approach the SIT spc;ctral gap 
for forward- and backward-propagating waves (Mantsyzov [1995]) in the limit of 
weak reflection. The allowed middle band corresponds to a polaritonic (collective 
atomic polarization) excitation in the periodic structure. It is different from single- 
photon hopping in a PBG via resonant dipole-dipole interactions (John and Quang 
[1995]). 



§ 5. Bright Solitons in RABR 

5.1. STANDING (QUIESCENT) SELF-LOCALIZED PULSES 

Stationary solutions of Eqs. (4.26a) and (4.26b) corresponding to bright solitons 
have been found by Kozhekin, Kurizki and Malomed [1998]. Such solutions for the 
symmetric-mode field E+ and polarization P are sought in the form 

E+ = e-'^^S{C), P = i e-'^^V{0 (5-1) 
with real V and <S. Substituting this into (4.26b), we eliminate V in favor of <S, 

and obtain an equation for SiC), 

S'' = ir,^-X^)S-2S ^^-f-''^^^-J\ (5.3) 

V (x - + ^ 

where the prime stands for d/d(. Equation (5.3) can be cast into the form of the 
Newton's equation of motion for a particle with the coordinate iS(C) moving in a 
potential U{S): 

S" = -U'{S), (5.4) 

where 



U{S) = --{rf - x^)S^ + 2(r? - x) ■ sign(x - 5)^{x-5Y+S\ (5.5) 
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The potential gives rise to bright solitons (Newell and Moloney [1992]), provided 
it has two symmetric minima (see Fig. 5). As follows from Eq. (5.5), the latter 
condition implies that the quadratic part of the potential is concave, i.e., |x| > r], 
and the second (asymptotically linear) part of the expression (5.5) is convex, so 
that X < S. Moreover, two minima separated by a local maximum in 5 = appear 
if U"{0) < 0. From this inequality it follows that bright solitons can appear in two 
frequency bands x, the lower band 

Xi < X < min{x2, S}, (5.6) 

and the upper band 

max{xi, r],5} <x< X2, (5.7) 
where the boundary frequencies Xi,2 are given by 



Xi,2 = (l/2) [d-r]TViv + Sf + 8\. (5.8) 

The lower band exists for all values ?7 > and 6, while the upper one only exists 
for 

S>V- 1/??, (5.9) 

which follows from the requirement X2 > V (see Eq. (5.7)). On comparing these 
expressions with the spectrum shown in Fig. 4, we conclude that part of the lower 
gap is always empty from solitons, while the upper gap is completely filled with 
stationary solitons in the weak-reflectivity case (5.9), and completely empty in the 
opposite limit. It is relevant to mention that a partly empty gap has also been found 
in a Bragg grating with second harmonic generation (Peschel, Peschel, Lederer and 
Malomed [1997]), see Ch. 2. 
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The bright soUton corresponds to the solution of the Newton equation (5.4) with 
the 'particle' sitting at time — oo on the local maximum 5 = 0, then swinging to 
one side and finally returning to 5 = at time +oo. Such solutions have been 
found in an implicit form by Kozhekin, Kurizki and Malomed [1998] : 



5(C) = 2|x- W)(i-^'(C))" 



with 



ICI = 



x-v 



(l-7e^)-i/2tan-i 



i-nl 



+(27eo)-^ In 



7^0 + ^nl - 7^2 



n 



and 



l(x + ^)(x-^)l 



(5.10) 



(5.11) 



(5.12) 



(note that TZq is positive under the above conditions (5.6)-(5.9)). It can be checked 
that this zero-velocity (ZV) gap soliton is always sm^/e-humped. Its amplitude can 
be found from Eq. (5.11), 



47eo/v^X + ??l- 



(5.13) 



The polarization amplitude V is determined by S via Eq. (5.2). 

To calculate the electric field in the antisymmetric E_ mode, we substitute 



into Eq. (4.21b) and obtain 



^" + (x' -ri^)A = 2V', 



(5.14) 



(5.15) 



which can be easily solved by the Fourier transform, once V{Q is known. An 
example of bright solitons is depicted in Fig. 6. Note that, depending on the 
parameters ry, 5 and x, the main part of the soliton energy can be carried either by 
the or the 5]_ mode. 

The most drastic difference of these new solitons from the well-known SIT pulses 
is that the area of the ZV soliton is not restricted to 27r, but, instead, may take an 
arbitrary value. As mentioned above, this basic new result shows that the Bragg 
reflector can enhance (by multiple reflections) the field coupling to the TLS, so as to 
make the pulse area effecUvely equivalent to 27r. In the limit of the small-amplitude 
and small-area solitons, T^q <g: 1, Eq. (5.11) can be easily inverted, the ZV soliton 
becoming a broad sech-like pulse: 



5«2|x-(5|7^osech 



X-V 



(5.16) 
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Figure 6: Zero-velocity (RABR) solitons |S+(C)|^ (a) S = Q, t] = 0.9, x = -0.901 
(divergent width and amplitude); (b) idem, but for x = 0.901 (divergent width and 
finite amplitude). 

In the opposite limit, 1 — TZq 0, i.e., for vanishingly small |x + ?/|, the soliton's 
amplitude (5.13) becomes very large, and further analysis reveals that, in this case, 

the soliton is characterized by a broad central part with a width ~ (l — TZq) 
(Fig. 6(a)). Another special limit is x — ^7 — * 0. It can be checked that in this limit, 
the amplitude (5.13) remains finite, but the soliton width diverges as |x — jyl^^^^ 
(Fig. 6(b)). Thus, although the ZV soliton has a single hump, its shape is, in 
general, strongly different from that of the traditional nonlinear-Schrodinger (NLS) 
sech pulse. 

5.1.1. Stability 

The stability of the ZV gap solitons was tested numerically, by means of direct 
simulations of the full system (4.26), the initial condition taken as the exact soliton 
with a small perturbation added to it. Simulations at randomly chosen values of 
the parameters, have invariably shown that the ZV GS are apparently stable. How- 
ever, the possibility of their dynamical and structural instability needs be further 
investigated, as has been done in the case of GS in a Kerr-nonlinear fiber with a 
grating by Barashenkov, Pelinovsky and Zemlyanaya [1998] and SchoUmann and 
Mayer [2000]. 

5.2. MOVING SOLITONS 

Although the system of Eqs. (4.26) is not explicitly Galilean- or Lorentz-invariant, 
translational invariance is expected on physical grounds. Hence, a full family of 
soliton solutions should have velocity as one of its parameters. This can be ex- 
plicitly demonstrated in the limit of the small-amplitude large-width solitons [cf. 
Eq. (5.16)]. We search for the corresponding solutions in the form Y,+ {(,t) — 
5(C,T)exp(-ixoT), P{(,t) = i'P(C, t) exp (-^xot) (cf. Eqs. (5.1)), where xo is 
the frequency corresponding to /c = on any of the three branches of the dispersion 
relation (4.58) (see Fig. 4), and the functions S((,t) and 'P{C,,t) are assumed to 
be slowly varying in comparison with exp (— ixoT). Under these assumptions, we 
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(b) 




Figure 7: Pulses obtained as a result of 'pushing' a zero-velocity RABR soliton 
(dashed lines): (a) push, characterized by the initial multiplier exp(— ipC^) after a 
sufficiently long evolution (r = 400) (solid lines). i5 = 0, ry = 4, x = —4.4, and 
p = 0.1. (b) idem, but for p = 0.5. 



arrive at the following asymptotic equation for S{C,t): 



Xo^V I 
(Xo-Sr 



S = [V'^xl + 2^]S. 



Since this equation is of the NLS form, it has the full two-parameter family of 
soliton solutions, including the moving ones (Newell and Moloney [1992]). 

In order to check the existence and stability of the moving solitons numerically, 
the following procedure has been used by Kozhekin, Kurizki and Malomed [1998]: 
Eqs. (4.26) were simulated for an initial configuration in the form of the ZV soliton 
multiplied by exp(zp^) with some wavenumber p, in order to 'push' the soliton. 
The results demonstrate that, at sufficiently small p, the 'push' indeed produces a 
moving stable soliton (Fig. 7(a)). However, if p is large enough, the multiplication 
by exp(zpC) turns out to be a more violent perturbation, splitting the initial pulse 
into two solitons, one quiescent and one moving (Fig. 7(b)). 

Another one-parameter subfamily of moving GS was found in the exact form of 
a phase-modulated 27r-soliton by Kozhekin and Kurizki [1995]: 

S+ = Ao exp [z («C - xr)] sech [/? (C - vt)] , (5.18) 

where x is the detuning from the gap center, Aq is the amplitude of the solitary 
pulse, P its width and v its group velocity. 

Substituting d^-P from Eq.(4.23a) into Eq.(4.21a), we may express P in terms of 
and the population inversion w. Then, upon eliminating P and using ansatz 
Eq.(5.18), we can integrate Eq. (4.23b) for the population inversion w, obtaining 



2{S - 77) cosh^ [/3(C - vt)] ' 



Using these explicit expressions for P and w in Eqs. (4.21a) and (4.23a), we reduce 
our system to a set of algebraic equations for the coefficients k, x that determine 
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the spatial and temporal phase modulation, and the pulse width (3 as functions of 
the velocity v. 

2{x-n/v)-{l-l/v^){5-ri) = 0, (5.20a) 

+20^v^{x- K/v) + 2{5-ri) = 0, (5.20b) 
(/3V-/32 + K2-x'+r?' + 2)-2(x + ^)(x-/«H = 0. (5.20c) 

The soliton amplitude is then found to satisfy |^o| = 2/3t), exactly as in the case 
of usual SIT (see Ch. 3). This implies, by means of Eq.(5.18), that the area under 
the 5]_|_ envelope is 27r. 

Let us consider the most illustrative case, when the atomic resonance is exactly 
at the center of the optical gap, 5 = 0. Then the solutions for the above parameters 
are 



2v ' 
1 

2Y 



(5.21a) 



X = (5.21b) 



8v^il -v^)-r^^{l + vY 
4(l-t;2)2 



/3V = |^o|74=^^^^ TTT • (5.21c) 



In the frame moving with the group velocity of the pulse, (' = (—vt, the temporal 

phase modulation will be (kv — x)t, which is foimd from Eq.(5.21) to be equal to 
—rjT. Since r] is the (dimensionless) 'bare' gap width (see Ch. 4, this means that 
the frequency is detuned in the moving frame exactly to the band-gap edge. The 
band-gap edge corresponds (by definition) to a standing wave, whence this result 
demonstrates that such a pulse is indeed a soliton, which does not disperse in its 
group-velocity frame. 

The allowed range of the solitary group velocities may be determined from 
Eq.(5.21c) through the condition /J^ > for a given 77. The same condition implies 
Ivl < ^max, where 



, _ 8.^(1-.^) 

'/max— ^l^y2^2 ' {O.ZZ) 



It follows from Eq.(5.22) that the condition for 2n SIT gap soliton (5.18) is \r]\ < 1, 
'ymax = 1 corresponding to ?; = 1/ This condition means that the cooperative 
absorption length ctq/uq should be shorter than the reflection (attenuation) length 
in the gap 4c/(aia;c'^o), i.e., that the incident light should be absorbed by the TLS 
before it is reflected by the Bragg structure. In addition, both these lengths should 
be much longer than the light wavelength for the weak-reflection and slow- varying 
approximation to be valid. 

From Eq. (4.21b) we find S_ = The envelopes of both waves (forward 

and backward) propagate in the same direction; therefore the group- velocity of the 
backward wave is in the direction opposite to its phase-velocity! This is analogous 
to climbing a descending escalator. 

Analogously to Kerr-nonlinear gap solitons (Ch. 2), the real part of the non- 
linear polarization ReP creates a traveling 'defect' in the periodic Bragg reflector 
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Figure 8: Dependence of the solitary pulse velocity (solid line) and amplitude 
(dashed line) in RABR on frequency detuning from the gap center for r/ = 0.7. 
The 'bare' gap edge marked by dotted line. 

structure which allows the propagation at band-gap frequencies. The real part of 
the nonlinear polarization is governed by the frequency detuning from the TLS res- 
onance. Exactly on resonance (which we here take to coincide with the gap center) 
X = 5 = 0, ReP = 0, and our solutions (5.21) yield imaginary values of the velocity 
V and modulation coefficient k. The forward field envelope then decays with the 
same exponent as in the absence of TLS in the structure. Because of this mecha- 
nism, SIT exists only on one side of the band-gap center, depending on whether the 
TLS are in the region of the higher or the lower linear refractive index. This result 
may be understood as the addition of a near-resonant non-linear 'refractive index' 
to the modulated index of refraction of the gap structure. When this addition com- 
pensates the linear modulation, soliton propagation becomes possible (sec Fig. 3). 
On the 'wrong' side of the band-gap center, soliton propagation is forbidden even in 
the allowed zone, because the nonlinear polarization then cannot compensate even 
for a very weak loss of the forward field due to reflection. 

The soliton amplitude and velocity dependence on frequency detuning from the 
gap center (which coincides with atomic resonance) are illustrated in Fig. 8. They 
demonstrate that forward soliton propagation is allowed well within the gap, for x 
satisfying (1 — ^/l — rf) /r/<x<(l + \/ 1 — ) /rj. In addition to frequency detun- 
ing from resonance, the near-resonant GS possesses another unique feature: spatial 
self-phase modulation kQ of both the forward and backward field components. 

5.3. NUMERICAL SIMULATIONS 

To check the stability of the analytical solution Eq.(5.18), as well as the possibility 
to launch a moving GS by the incident light field, numerical simulations of Eqs.(4.6) 
were performed by Kozhekin and Kurizki [1995]. As the launching condition, the 
incident wave was taken in the form £f — Acxp [ix(t — to)]/ cosh [/3(t — to)/''o] 
without a backward wave {£b = 0) at the boundary of the sample z = 0. By 
varying the detuning x and amplitude A we investigate the field evolution inside the 
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Figure 9: Numerical simulations of the intensities of (a) 'forward' and (b) 'back- 
ward' waves in the RABR gap, when Eqs.(5.21), (5.22) are obeyed (tj = 0.7, group 
velocity f ~ 0.3) 



structure. When these parameters are close to those allowed by Eqs.(5.21), (5.22), 
we observe the formation and lossless propagation of both forward and backward 
soliton-like pulses with amplitude ratios predicted by our solutions (Fig. 9). By 
contrast, exponential decay of the forward pulse in the gap is numerically obtained 
in the absence of TLS (Fig. 10). 

The analysis surveyed in Sees. 5.1-5.3 strongly suggests, but does not rigor- 
ously prove, that the solution subfamily (5.18) belongs to a far more general two- 
parameter family, whose other particular representatives are the exact ZV solitons 
(5.11) and the approximate small-amplitude solitons determined by Eq. (5.17). 



5.4. COLLISIONS BETWEEN GAP SOLITONS 

An issue of obvious interest is that of collisions between GSs moving at different 

velocities in RABR. In the asymptotic small-amplitude limit reducing to the NLS 
equation (5.17), the collision must be elastic. To get a more general insight, we 
simulated collisions between two solitons given by (5.18). The conclusion is that 
the collision is always inelastic, directly attesting to the nonintegrability of the 
model. Typical results are displayed in Fig. 11, which demonstrates that the 
inelasticity may be strong, depending on the parameters. 
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Figure 11: Typical example of inelastic collisions between the RABR solitons (5.18) 
at (5 = and rj = 0.5, with the velocities (normalized to c) ui — 0.6, U2 — —0.75. 
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§ 6. Dark solitons in RABR 

6.1. EXISTENCE CONDITIONS AND THE FORM OF THE SOLITON 

Dark solitons (DSs) in RABR have been studied by Opatrny, Malomed and Kurizki 
[1999]. They are obtained similarly to the bright ones, by solving Eq. (5.4) with 
the potential (5.5). The potential will give rise to DS's provided that it has two 
symmetric maxima (see Fig. 5). In this case the quadratic part of the potential is 
convex, i.e., |x| < and the second (asymptotically linear) part of the expression 
(5.5) is concave, so that x > ^- From these two inequalities, a simple necessary 
restriction on the model's parameters follows, 

S<r]. (6.1) 



The condition for the existence of the symmetric maxima determines the following 
frequency interval x (recall that r] is defined to be positive): 

max{(5, -77} <x< min{x2, v}, (6-2) 

Making use of (5.8), one can easily check that, once the condition (6.1) is satisfied, 
the DS-supporting band (6.2) always exists. The DS frequency range defined by 
Eq. (6.2) is marked by shading (to the right from zero) in Fig. 4. The maxima of 
the potential are located at the points 



Sm = ±V4(r? + x)-2-(x-^)^ (6.3) 
which correspond to the polarization values 

Vm = TV^-imix-Snx + V?- (6.4) 

Integrating Eq. (5.4) by means of energy conservation in the formal mechanical 
problem, we obtain S{Q in an implicit form. 



c 



±j,r I (6.5) 



with 



a='^{rf-x^),l3 = 2{r]-x),l = X-5. (6.6) 

The solution (6.5) corresponds to a trajectory beginning at the potential maximum 
±5m at 'time' C = —00 and arriving at the other maximum, ^Sm at 'time' C, = 00 
(see Fig. 5). In terms of the S+ mode of the electric field, this is exactly a quiescent 
(zero- velocity) DS with the background cw amplitude Sm- 

The integral (6.5) can be formally expressed in terms of incomplete elliptic inte- 
grals, but, practically, it is more helpful to evaluate it numerically. As in Ch. 5, 
the polarization amplitude V is determined by <S via Eq. (5.2) and the amplitude 
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Figure 12: A typical example of a dark soliton, presented in terms of the variables 
<S, V and A. The parameters are r\ = 0.6, (5 = —2, and x = 0.25 

of the S_ mode is obtained by solving Eq. (5.15). An example of the amplitude S 
for DS in the mode, together with the corresponding quantities V and A^ are 
plotted in Fig. 12. 

The energy density of the field mode always has the shape of a hole in the 

background (see Fig. 13). The energy density of has a hump, which is the 

counterpart of the hole in the mode. The net electromagnetic energy density 
may have either a hole (which never drops to zero) or a hump, depending on the 
system parameters ?7, 5, and the soliton frequency x- 

6.2. THE BACKGROUND STABILITY 

An obvious necessary condition for the stability of DS is the stability of its cw 
background. To tackle this problem, we use Eq. (4.26b) to eliminate S+ in favor 
of P, 

E+ = -(P.-^5P)(l-|P|2)-^/^ (6.7) 

and insert it into Eq. (4.26a). The resulting equation for P is linearized around 
the stationary value Vm (see Eq. (6.4)), substituting 

P = Vue-'^^ [1 + a(C, r) + i b{C, r)] , (6.8) 

where a and b are small real perturbations. We can look for its general solution in 
the form (cf. Eq. (4.56c)), 

a(C, r) = ooe'(««-"^) , 6(C, r) = 6oe'(««-"^) , (6.9) 

which leads to a dispersion relation for Cl and k, that consists of two parts: 

-(x-5)(x+^?)'^^' - 2 [x{x-s){x+vr + 2] 

+(??'-x') [4- (x-(5)'(r?+x)'] +4k' = 0, (6.10) 
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Figure 13: Field energy densities as a function of the coordinate ( inside the dark 
soliton shown in Fig. 12: (dashed hne), (dash-dotted hne), and their 

sum (continuous line). 

and 

-ft^ - (3x-<5)n' + [{v-x){2x+v-S)+K^] n + (x-sy = o. (e.ii) 

As checked numerically by Opatrny, Malomcd and Kurizki [1999] for many values 
of r], 6 and x that support DS's according to the results obtained in the previous 
section, all the roots of (6.10) and (6.11) are real for any real k. This implies the 
stability of the background for these values of 77, S and x- 

Equations (6.10) and (6.11) represent new dispersion relations, which are valid 
under the condition of strong background field Sm and replace the zero-field dis- 
persion relations of Eq. (4.57). This kind of optical bistability can be compared 
to the distributed feedback bistability with Kerr nonlinearity studied by Winful, 
Marburger and Garmire [1979]. 

6.3. DIRECT NUMERICAL STABILITY TESTS 

Even though there is no evidence of background instability, it is necessary to simu- 
late the full system of the partial differential equations, in order to directly test the 
DS stability. Equations (4.26a) and (4.26b) have been integrated numerically by 
Opatrny, Malomed and Kurizki [1999], with initial conditions differing from the ex- 
act DS solution by a small perturbation added to it. The results strongly depend on 
the parameters r], 6 and x- for some values, an explosion of the initial perturbation 
occurs, leading to a completely irregular pattern, whereas for others the DS shape 
remained virtually undisturbed. The dependence of the stability on the parameters 
T] and X with fixed 6 is shown in Fig. 14. The darkest area of the DS parameter 
region corresponds to the stable regime where no instability has occurred during 
the entire simulation time (typically, r ~ 500). In the rest of the parameter region, 
DS's are unstable: in the lightest part of the DS region, the instability develops 
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Figure 14: Parameter regions (ry vs. x) for dark (Eq. (6.2)) and bright (Eqs. (5.6) 
and (5.7)) solitons a,t S — —2. The boundaries xi and X2 are given by Eq. (5.8). 
In the DS region, the darkest area corresponds to stable behavior whereas in the 
remaining part the (numerical) solutions are unstable: in the lightest area (of DS's) 
instability develops very quickly, while in the intermediate area the DS exists for a 
much longer time before the onset of the instability. 

very quickly (at r < 50), whereas in the intermediate part, the instability builds up 
relatively slowly. As can be seen, the unstable behavior occurs closer to the bound- 
aries of the existence region, x = ^ V and x = whereas along the boundary x 
= X2 (corresponding to the DS-supporting background which degenerates into the 
trivial zero solution), DS's are stable. 

6.4. COEXISTENCE OF THE DARK AND BRIGHT SOLUTIONS 

A very interesting question is whether the system can support bright and dark 
solitons at the same values of the parameters. As mentioned above, DS's always 
exist in the frequency interval (6.2), once the inequality (6.1) is satisfied. On the 
other hand, bright solitons arc found in two frequency bands x, given in Eqs. (5.6) 
and (5.7). From the discussion in Ch. 5 it follows that the DS frequency band 
always coexists with one or two bands supporting the bright solitons. The special 
case when there arc two bright-soliton bands coexisting with the DS band is singled 
out by the condition 



One can readily check that the coexisting frequency bands supporting bright and 
dark solitons never overlap, i.e., quite naturally, the bright and dark solitons cannot 
have the same frequency. 

6.5. MOVING DARK SOLITONS 

Thus far we have considered only the quiescent DS's. A challenging question is 
whether they also have their moving counterparts. Adding the velocity parameter 
to the exact DS solution is not trivial, as the underlying equations (4.26a) and 
(4.26b) have no Galilean or Lorentzian invariance. The physical reason for this is 



T] < S < ri. 



(6.12) 
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Figure 15: Moving dark soliton: the values of the parameters are the same as in 
Fig. 12, the background phase-jump parameter (f) (see Eqs. (6.13) and (6.14)) is 
<f) = — 7r/4. Dashed line: r = 0, continuous line: r = 600. 



the existence of the special (laboratory) reference frame, in which the Bragg grating 
is at rest. 

In principle it is possible, in analogy to the stationary solutions and Eq. (5.1), sub- 
stitute functions of the argument (( — vt) into the set (4.26a) and (4.26b) to obtain 
an ordinary differential equation. However, this would be a complicated complex 
nonlinear equation of the third order, containing all the lower-order derivatives, so 
that we would not be able to take advantage of the Newton-like structure, as in 
Ch. 5 and 6. Though it is possible to solve such an equation numerically, it is more 
suitable to deal with the original set of partial differential equations, in order to 
better understand the nature of the evolution. 

In contrast to the case of bright solitons where the moving solutions can be found 
by multiplying, in the initial conditions, the quiescent DS by a factor proportional 
to exp(iKC) (see Kozhekin, Kurizki and Malomed [1998] and Ch. 5), it has proven 
possible to generate stable moving DS's from the quiescent ones in a different way 
(Opatrny, Malomed and Kurizki [1999]). To this end, recall that a DS corresponds 
to a transition between two different values of the background cw field. The back- 
ground field takes, generally, complex values (note the real values in the expressions 
(6.3) and (6.4) are only our choices adopted above for convcnicnc;e). The quiescent 
DS corresponds to a transition between two background values with phases differ- 
ing by TT. A principal difference of the DS's in the present model from those in 
the NLS equation (Kivsliar and Luther-Davies [1998]) is that here a moving DS is 
generated by introducing a phase jump ^ tt across the DS. 

Thus, one can take the initial condition for the system of equations (4.26a) and 
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(4.26b) as 

S+(C,0) = cos[ ^ ]Sg{0+ ism {^]Sm, (6.13) 



P = cos(^^^Pg(C)+isin(^0pM, (6.14) 

where Sq and Vq are the (real) functions corresponding to the quiescent DS, Sm and 

Vm arc given by Eqs. (6.3), (6.4) and </> is the deviation from tt of the background 
phase jump across DS. A typical result obtained by means of this modification of the 
initial state is displayed in Fig. 15: the DS moves at a velocity that is proportional 
to 4>. The rcsiilting form of the moving DS is slightly different from that of the 
quiescent soliton. The moving DS appears to be stable over the entire simulation 
time. 

§ 7. Light bullets (spatiotemporal solitons) 

A promising direction is the study of solitons in resonantly absorbing multi-dimensional 
(2D and 3D) media, in which the quasi-lD periodic structures can be realized as 
thin homogeneous layers set perpendicular to the direction of light propagation. 
In such media, spatiotem,poral solitons, i.e., those localized in all dimensions, both 
transverse (spatial proper) and longitudinal (effectively, temporal), may exist. Spa- 
tiotemporal optical solitons or 'light bullets' (LBs) in various nonlinear media are 
surveyed in Ch. 1. 

Here we are concerned with LBs in RABRs that consist of thin TLS layers em- 
bedded in a 2D- or 3D- periodic dielectric medium. We will follow a recent analysis 

by Blaaubocr, Kurizki and Malomed [2000], which extends an earlier prediction of 
stable LBs in uniform 2D and 3D SIT media (Blaauboer, Malomed and Kurizki 
[2000]). 

We start by considering a 2D SIT medium with a refractive index n{z, x) peri- 
odically modulated in the propagation direction z, which represents the quasi-one- 
dimensional Bragg grating. Light propagation in the medium is described by the 
lossless Maxwell-Bloch equations (Newell and Moloney [1992]): 

dP 

— = 0, (7.1b) 
^ + h£*P + P*£) = 0. (7.1c) 

OT 2 

Here (as in Sec. 4.1) £ and P are the slowly varying amplitudes of the electric field 
and medium's polarization, w is the population inversion, ( and x are longitudinal 
and transverse coordinates (measured in units of the resonant-absorption length), 
and T is time (measured in units of the input pulse duration). The Fresnel num- 
ber, which governs the transverse diffraction in the 2D and 3D propagation, was 
incorporated into x, and the detuning of the carrier frequency luq from the central 
atomic-resonance frequency was absorbed into £ and P. To neglect the polariza- 
tion dephasing and inversion decay, we assume pulse durations that are short on 
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the time scale of the relaxation processes. Equations (7.1) are then compatible with 
the local constraint \P\'^ + ~ 1, which represents the so-called Bloch- vector con- 
servation. In a ID case, i.e., in the absence of the a;-dependence and for n{z, x) = 1, 
Eq. (7.1a) reduces to the sine-Gordon (SG) equation, which has a commonly known 
soliton solution, see Eq.(3.11) and Ch.3. 

To search for LBs in a 2D medium subject to a resonant periodic longitudinal 
modulation, one may assume a periodic modulation of the refractive index as per 
Eq.(4.1). The RABR is then constructed by placing very thin layers (much thinner 
than 1/kc) of two- level atoms, whose resonance frequency is close to Wc, at maxima 
of this modulated refractive index. 

The objective is to consider the propagation of an electromagnetic wave with a 
frequency close to Wc through a 2D RABR. Due to the Bragg reflections, the electric 
field £ gets decomposed into forward- and backward-propagating components £p 
and £b, which satisfy equations that are a straightforward generalization of the ID 
equations derived by Kozhekin and Kurizki [1995], Kozhekin, Kurizki and Malomcd 
[1998], and Opatrny, Malomed and Kurizki [1999] (see also Eqs. (4.21) and (4.23) 
in this review): 

(7.2a) 

0, 

(7.2b) 

0, 

0, (7.2c) 
0. (7.2d) 

Here S± is defined by Eq. (4.9), and 77 is a ratio of the resonant-absorption length 
in the two-level medium to the Bragg-reflection length, which was defined above by 
Eq. (4.13). 

To construct an analytical approximation to the LB solutions, the starting point 
adopted by Blaauboer, Kurizki and Malomed [2000] is a subfamily of the exact ID 
soliton solutions to Eqs. (7.2), which was found by Kozhekin and Kurizki [1995] (see 
also Sec. 5.2 in this review) and is given by Eqs. (5. 18) and (5.19). These solutions 
were taken with parameters satisfying Eqs. (5.20): 

Ao = 2^, P=^l, - = -/f^> (7-3a) 
K = — \/p^^rf, and x = (5. (7.3b) 

These solutions were chosen as a pattern to construct an approximate solution 
for LBs because the shape of the fields and S_ in the solutions is similar to that 
of the SG soliton in the ID uniform SIT medium (see Sec. 3.1). Inspired by this 
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analogy and by the fact that there exist LBs in the uniform 2D SIT medium which 
reduce to the SG solitons in the ID limit (Blaauboer, Malomed and Kurizki [2000]), 
one can search for an approximate LB solution to the 2D equations (7.2), which 
also reduces to the exact soliton in ID. To this end, the following approximation 
was assumed: 



with ei(T, C) ^ /3(C - vt) + 60 + Cx, 62 (r, C) = /?(C - vt) + Go - Cx, the phase 
v and coefficients Qq and C being real constants, while the other parameters are 
defined by Eqs.(7.3). 

The ansatz (7.4) satisfies Eqs. (7.2a) and (7.2b) exactly, while Eqs. (7. 2d) are 
satisfied to order y/5/r] — IC^, which requires that y^S/rj — IC^ <C 1. The ansatz 
applies for arbitrary r], admitting both weak {rj <^ 1) and strong {r] > 1) reflectivi- 
ties of the Bragg grating, provided that the detuning remains small with respect to 
the gap frequency. Comparison with numerical simulations of Eqs. (7.2), using Eq. 
(7.4) as an initial configuration (a finite-difference method, with Fourier transform 
scheme, described by Drummond [1983], was used), tests this analytical approxi- 
mation and shows that it is indeed fairly close to a numerically exact solution; in 
particular, the shape of the bullet remains within 98% of its originally presumed 
shape after having propagated a large distance, as is shown in Fig. 16. 

Three-dimensional LB solutions with axial symmetry have also been constructed 
in an approximate analytical form and succesfully tested in direct simulations, 
following a similar approach (Blaauboer, Kurizki and Malomed [2000]). Generally, 
they are not drastically different from their 2D counterparts described above. 

A challenging problem which remains to be considered is the construction of 
spinning light bullets in the 3D case (doughnut-shaped solitons, with a hole in 
the center, carrying an intrinsic angular momentum). Recently, spinning bullets 
were found by means of a sophisticated version of the variational approximation 
in a simpler 3D model, viz., the nonlinear Schrodinger equation with self- focusing 
cubic and self-defocusing quintic nonlinearities, by Desyatnikov, Maimistov and 
Malomed [2000]. Further direct simulations have demonstrated that these spinning 
bullets (unlike their zero-spin counterparts) are always subject to an azimuthal 
instability, that eventually splits them into a few moving zero-spin solitons, although 
the instability can sometimes be very weak (Mihalache, Mazilu, Crasovan, Malomed 
and Ledercr [2000]). At present, it is not known whether spinning LBs can be 
completely stable in any 3D model. 

It is relevant to stress that two- (and three-) dimensional LB solutions of the 
variable-separated form, S+ ^ S_ ^ .f('''X) ' 3(2;), do not exist in the RABR 
model. Indeed, the substitution of this into Eqs. (7.2a) and (7.2b) yields only a 
plane- wave solution of the form S± exp {iAr + iBx), with constant A and B. 



S_ = 



P = 
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Figure 16: The forward-propagating electric field of the two-dimensional 'light bul- 
let' in the Bragg reflector, \£f\, vs. time t and transverse coordinate x, after having 
propagated the distance z = 1000. The parameters are 77 — 0.1, 6 — 0.2, C = 0.1 
and 80 = —1000. The field is scaled by the constant h/ATQiinQ. 

§ 8. Experimental prospects and conclusions 

This review has focused on properties of solitons in RABR, combining a periodic 
refractive-index (Bragg) grating and a periodic set of thin active layers (consisting of 
two-level systems resonantly interacting with the field). It has been demonstrated 
that the RABR supports a vast family of bright gap solitons, whose properties 
differ substantially from their counterparts in periodic structures with either cubic 
or quadratic off-resonant nonlinearies reviewed in Ch. 2. 

The same RABR can support, depending on the initial conditions, either dark 
or bright stable solitons, without any changes of the system parameters, which is 
a unique feature for nonlinear optical media (Ch. 5, 6). Zero-velocity dark soli- 
tons can be found in an analytical form, as well as traveling dark solitons with a 
constant phase difference (7^ tt) of the background amplitudes across the soliton 
(Ch. 6). The latter property is a major difference with respect to dark solitons of 
the NLS equation, whose motion is supported by giving the background a nonzero 
wavenumber. 

Depending on the values of the parameters, the frequency band of the quiescent 
dark solitons coexists with one or two bands of the stable bright ones, without an 
overlap. Direct numerical simulations demonstrate that some dark-soliton solutions 
are stable against arbitrary small perturbations, whereas others are unstable when 
they are close to the "dangerous" boundaries of their existence domain. 

A multidimensional version of the RABR model, corresponding to a periodic 
set of thin active layers placed at the maxima of the refractive index, which is 
modulated along the propagation direction of light has been considered too. It has 
been found to support stable propagation of spatiotemporal solitons in the form of 
two- and three-dimensional 'light bullets' (LBs). 

The best prospect of realizing a RABR which is adequate for observing the soli- 
tons and light bullets discussed in Ch. 5 - Ch. 7 is to use thin layers of rare-earth 
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ions (Greiner, Boggs, Loftus, Wang and Mossbcrg [1999]) embedded in a spatially- 
periodic semiconductor structure (Khitrova, Gibbs, Jahnke, Kira and Koch [1999]). 
The two-level atoms in the layers should be rare-earth-ions with the density of 
10^^ — 10^® cm~'^, and large transition dipole moments. The parameter rj can vary 
from to 100 and the detuning is ~ 10^^ — 10^'^ s~^. 

Cryogenic conditions in such structures can strongly extend the dephasing time 
T2 and thus the soliton's or LB's lifetime, well into the /isec range (Greiner, Boggs, 
Loftus, Wang and Mossberg [1999]), which would greatly facilitate the experiment. 
The construction of suitable structures constitutes a feasible experimental challenge. 

In a RABR with the transverse size of 10 /im, LBs can be envisaged to be localized 
on the time and transverse-length scales, respectively, ~ 10~^^ s and 1 /Ltm. The 
incident pulse has uniform transverse intensity and the transverse diffraction is 
strong enough. One needs rf^/Zabs-^o < 1, where Zabs. -^o and d are the resonant- 
absorption length, carrier wavelength, and the pulse diameter, respectively (Slusher 
[1974]). For l^hs ~ 10~^ m and Aq ^ 10~^ m, one thus requires d < 10~^ m, which 
implies that the transverse medium size must be a few /im. 

Effects of TLS dephasing and deexcitation in RABRs can be studied by substi- 
tuting the values —i6—r2TQ for the frequency term —iS in Eqs. (4.47) - (4.49) and 
the loss terms -riTo(wo + 1) in Eq. (4.50), -Tiwi in Eq. (4.51) and -ri(ui2 - 2) 
in Eq. (4.52). We have checked that these modifications do not influence the qual- 
itative behavior of the solutions on the time scale ttq < l/ri,2. 

Let us now discuss the experimental conditions for the realization of RABR 
solitons using quantum wells embedded in a semiconductor structure with peri- 
odically alternating linear index of refraction (Khitrova, Gibbs, Jahnke, Kira and 
Koch [1999]). We can assume the following values: the average refraction index is 
no ~ 3.6, the wavelength (in the medium) A « 232 nm, which corresponds to the 
angular frequency cJc ~ 2.26x10^^ s~^. Excitons in quantum wells can, under cer- 
tain conditions (such as low densities and proximity of the operating frequency to 
an excitonic resonance, see Khitrova, Gibbs, Jahnke, Kira and Koch [1999]) may be 
regarded as effective two-level systems (TLS's). We consider their surface density 
to be « 10-^° - 10^^ cm"^, which corresponds to a bulk density po ~ 10"^^ - 10^^ 
cm~^. If we assume that the excitons are formed by electrons and holes displaced 
by « 1 - 10 nm, then the characteristic absorption time Tq defined in Eq. (3.8) is 
To « 10""'^^ - lO"'^^ s, and the corresponding absorption length is cro/no w 10 - 
100 /xm. The dephasing time for excitons discussed by Khitrova, Gibbs, Jahnke, 
Kira and Koch [1999] is l/r2 « 10~^^ s, which seems to be the chief limitation 
of the soliton lifetime for this system. The structures shown in Figs. 6, 12, 13, 
and 15, occupying regions of approximately 100 absorption lengths would require 
a device of the total width of approximately 1 mm to 1 cm, which corresponds to 
PS 10'^ to 10^ unit cells. The modulation of the refraction index can be as high as 
ai « 0.3, so that the parameter rj (see Eq. (4.13)) can vary from to 100. The 
unit of the dimensionless detuning 6 would represent a 10~^ - 10~^ fraction of the 
carrier frequency. The intensities of the applied laser field corresponding to E± w 1 
are then of the order 10^ - 10'' W/cm^. 

In the work by Khitrova, Gibbs, Jahnke, Kira and Koch [1999], the width of the 
active layers (quantum wells) is considered to be 5 - 20 nm, which corresponds to 
the parameter 7^ (see Eq. (4.29)) in the range 10^"^ - 2 x 10~^. In the simulations 
discussed by Opatrny, Malomed and Kurizki [1999], taking the largest of these 
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values and the parameters as in Fig. 12, i.e., rj = 0.6, S = —2, and x — 0.25, we 
have observed the time evolution of the system (4.36a) (4.52). As the initial 
condition, both the DS solution corresponding to zero width of the active layers 
and the DS solution including the finite width correction, have been taken. In both 
cases, the evolution was quite regular over the observed time r 50, and the zero- 
width solution (with the quantities S,'P,A as given by Eqs. (6.5), (5.2) and (5.15)) 
started to change after t « 10. 

We can now sum up the discussion of experimental perspectives for the realization 
of RABR solitons: a) The prospects appear to be good for gratings incorporating 
thin layers of rare-earth ions under cryogenic conditions, b) The realization of these 
solitons in excitonic superlattices would require much longer dephasing times than 
those currently achievable in such structures. 
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